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O ■ ABSTRACT 



We have made timing observations of binary pulsar PSR B1534+12 with radio 
telescopes at Arecibo, Green Bank, and Jodrell Bank. By combining our new 
Q ■ observations with data collected up to seven years earlier, we obtain a significantly 

CS| ■ improved solution for the astrometric, spin, and orbital parameters of the system. 

For the first time in any binary pulsar system, no fewer than five relativistic 
or "post-Keplerian" orbital parameters are measurable with useful accuracies 
in a theory-independent way. We find the orbital period of the system to be 



decreasing at a rate close to that expected from gravitational radiation damping, 
according to general relativity, although the precision of this test is limited to 
about 15% by the otherwise poorly known distance to the pulsar. The remaining 
Q\ ■ post-Keplerian parameters are all consistent with one another and all but one 

of them have fractional accuracies better than 1%. By assuming that general 



relativity is the correct theory of gravity, at least to the accuracy demanded by 
this experiment, we find the masses of the pulsar and companion star each to be 
1.339 ± O.OO3M and the system's distance to be d = 1.1 ± 0.2 kpc, marginally 
larger than the d ~ 0.7 kpc estimated from the dispersion measure. The increased 



distance reduces estimates of the projected rate of coalescence of double neutron- 
star systems in the universe, a quantity of considerable interest for experiments 
with terrestrial gravitational wave detectors such as LIGO. 

Subject headings: pulsars: individual (PSR B1534+12) — gravitation — binaries: 
close — stars: distances 
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1. Introduction 



Pulsars in double-neutron-star binaries provide the best known laboratories for ex- 
perimental tests of gravity in the radiative and strong field regimes. Timing analysis of 
pulsar signals allows measurement of five Keplerian orbital elements as well as a number 
of post-Keplerian (PK) orbital parameters. The PK parameters can be analyzed using a 
theory-independent procedure in which the masses of the two stars are the only dynamically 
important a priori unknowns ( Pamour fc Taylor 1992 ). Each of the PK parameters depends 



on the masses in a different way; consequently, if any two of them are measured, the relevant 
parameters of the two-body system are fully determined within any gravitational theory. If 
three or more PK parameters can be measured, the overdetermined system can be used to 
test the gravitational theory itself. Measured values of the PK parameters Co (rate of advance 
of periastron), 7 (time dilation and gravitational redshift) and P& (orbital period derivative) 
for binary pulsar PSR B1913+16 have been found to be in excellent agreement with the 
predictions of general relativity ( Taylor fc Weisberg 1989] , pamour fc Taylor 1991| , Taylor 



|1994|) . In particular, the observed value of Pb confirms that the system is losing energy in the 
form of quadrupolar gravitational waves at the rate expected according to general relativity 
(GR). 

An opportunity to repeat and extend this test of relativistic gravitation was provided by 
discovery of a second suitable binary pulsar, PSR B1534+12, in a 10.1 hour orbit ( |Wolszczan] 
[19911) . PSR B1534+12 is significantly brighter than PSR B1913+16, and its pulse has a 
narrow peak, allowing more precise timing measurements. Because the orbit is nearly edge- 
on as viewed from the Earth, the PK parameters r and s (the "range" and "shape" of the 
Shapiro time delay) are more easily measured. In addition, both Co and 7 are measurable, 
as for PSR B1913+16, because of their contribution to secularly accumulating effects in 
the data. The resulting overdetermination of the orbit has already led to a non-radiative 
test of gravitation theory in the strong-field regime, complementing the Cj-j-Pi, test for 
PSR B1913+16 ( llaylor et oTTggg ). 



Previously published timing measurements of PSR B1534+12 were all made with the 



Arecibo 305 m telescope. We have extended this sequence of observations flArzoumanian 



1995|) through early 1994, when the telescope went out of normal service for a major up- 
grading; since then we have used the 43 m telescope at Green Bank, West Virginia, and the 
76 m Lovell telescope at Jodrell Bank, England, to acquire additional data. Owing to its 
slower and less eccentric orbit, the expected rate of orbital period decay for PSR B1534+12 
is more than an order of magnitude less than for PSR B1913+16. However the timing data 
are considerably more precise, and we find that fitting our full sequence of pulse times of 
arrival (TOAs) to the standard relativistic pulsar timing model now yields a measurement 
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of Pf, in the solar-system barycentric rest frame with uncertainty approximately 7% of the 
expected GR value. The result is a convincing second test of the dissipative coupling be- 
tween accelerating masses and a gravitational radiation field that carries energy and angular 
momentum away from the system. Like the first such test ( |l'aylor fc Weisberg 1989 , 1'aylor 



1994|) these results are in good accord with general relativity theory. The precision of the 
new result is limited by uncertainty in the necessary kinematic corrections, which depend on 
the poorly known pulsar distance. Interestingly, we can invert the test, assuming that GR 
is the correct theory of gravity, and calculate that the distance of PSR B1534+12 must be 
about 1.6 times that estimated from the dispersion measure ( [Taylor fc Gordes 1993|) . 



2. Observations 

Observations at the 305 m Arecibo telescope were made primarily with the Princeton 
Mark III system ( [Stinebring et al. 1992| ), using dual-polarization 2 x 32-channel filterbanks 
with 0.25 MHz channels at 430 MHz and 1.25 MHz channels at 1400 MHz. The signals in 
each channel were square-law detected, and the two polarizations summed in hardware and 
then folded at the topocentric pulsar period. Integration times were typically 3 minutes 
at 430 MHz and 5 minutes at 1400 MHz. The resulting total-intensity profiles have 512 
(430 MHz) or 1024 (1400 MHz) bins across the 37.9 ms period. The effective time resolution 
is dominated by the dispersion smearing in a single channel, approximately 304 [is at 430 MHz 
and 44 /is at 1400 MHz; quadrature addition of the post-detection time constants and the 
boxcar-averaging bin sizes used for these observations yields the total resolutions quoted in 
Table |l|. The Arecibo observations extend from August 1990 through March 1994. 

We observed PSR B1534+12 with the National Radio Astronomy Observatory 43 m 
telescope at Green Bank West Virginia, between March 1994 and May 1997. Observing 
sessions were spaced typically at two month intervals, and we used frequencies of 575 and 
800 MHz. A digital Fourier transform spectrometer (the "Spectral Processor") analyzed 
signals in each of two polarizations into 512 spectral channels across a 40 MHz passband. 
The spectra were folded synchronously at the predicted topocentric pulse period and averaged 
for three minutes. The resulting pulse profiles had 128 phase bins, or a resolution of 296 /is. 
The data were de-dispersed after detection and opposite polarizations summed to produce 
a single total-intensity pulse profile for each integration. 

We observed PSR B1534+12 with the 76 m Lovell Telescope between January and July 
1997, using the Princeton Mark IV data acquisition system at 610 MHz ([Stairs, Taylor, <£] 
Thorsett 1998 ). Exploratory observations were made in January and February, followed 



by a concentrated campaign from 16 June to 14 July during which time the pulsar was 
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observed for up to 3 hours, during scintillation maxima, nearly every day. For each sense 
of circular polarization, a 5 MHz bandpass was mixed to baseband using local oscillators 
in phase quadrature. The four resulting signals were low-pass filtered at 2.35 MHz (filters 
60 dB down at 2.5 MHz), sampled at 5 MHz, quantized to 4 bits, and written to a large 
disk array and later to magnetic tapes. Upon playback, the undetected signal voltages were 
dedispersed using the phase-coherent technique described by Hankins and Rickett (1975). 
After amplitude calibrations are applied, self- and cross-products of the right- and left- 
handed complex voltages yield the Stokes parameters of the incoming signal. These products 
were folded at the topocentric pulsar period using 2048 phase bins, and pulse TOAs were 
determined from the resulting total-intensity profiles using 190 s integrations. A summary 
of the more important parameters and statistics of all four observing systems is presented 
in Table [TJ 

We used the same TOA-fitting procedure for all data sets. Each observed profile was 
fitted to a standard template, using a least-squares method in the Fourier transform domain 
to measure its time offset ( [Taylor 1992|) . The offset was added to the time of the first sample 



of a period near the middle of the integration, thereby yielding an effective pulse arrival 
time. A different standard template was used for each observing system and frequency; 
they were made by averaging the available profiles over several hours or more. Uncertainties 
in the TOAs were estimated from the least squares procedure, and also from the observed 
scatter of the TOAs within 30 minutes of each one. Each observatory's local time standard 
was corrected retroactively to the UTC timescale, using data from the Global Positioning 
System (GPS) satellites. 



3. Data Analysis 

3.1. The Timing Model 

A pulse received on Earth at topocentric time t is emitted at a time in the comoving 
pulsar frame given by 

T = t-t + Ac- D/f + A m + A EQ - A SQ - A fi - A £ - A 5 . (1) 

Here to is a reference epoch and Ac is the offset between the observatory master clock 
and the reference standard of terrestrial time. The dispersive delay is D/f 2 , where D = 
DM/2.41 x 10 -4 , with dispersion measure DM in cm _3 pc, radio frequency / in MHz, and 
the delay in seconds. Finally, A# , A^ , and A^ are propagation delays and relativistic 
time adjustments for effects within the solar system, and A#, A E and As are similar terms 
accounting for phenomena within the pulsar's orbit ( pamour fc Deruelle 198"6| , [Taylor & 
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Weisberg 1989, Damour fc Taylor 1992| ). (We ignore the uninteresting constant and uniform 
rate of change of the overall propagation delay.) The orbital A terms are defined by: 

Ar = xsmu>(cosu — e) + x(l — e 2 ) 1 / 2 cosujsmu, (2) 
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At a given time t, then, the propagation delay across the pulsar orbit is calculated by a 
model which incorporates ten parameters implicitly defined in the above equations: five 
Keplerian parameters (x, u>, T , Pf,, e) and five PK parameters (u>, Pb, 7, r, s). These 
quantities, in conjunction with a simple time polynomial to model the spin of the pulsar and 
with astrometric parameters to model the propagation of the signal across the solar system, 
constitute the free parameters to be fit in the theory-independent timing model. 

In a particular theory of gravity, the five PK parameters can be written as functions 
of the pulsar and companion star masses, mi and 7712, and the well-determined Keplerian 
parameters. Of particular interest is general relativity, where the equations are as follows 
(see P amour fc Deruelle 1986| , Taylor fc Weisberg 1989| , |Damour fc Taylor 1992| ): 
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(12) 



Here the masses mi, m2, and M = mi + m 2 are expressed in solar units, and we use the 
additional shorthand notations s = sin i and T = GM & /c 3 = 4.925490947 /xs, where i is 
the angle between the orbital angular momentum and the line of sight, G the Newtonian 
constant of gravity, and c the speed of light. 
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3.2. Arrival Time Analysis 



We used the standard tempo analysis software flTaylor fc Weisberg 1989| see also 



Internet location http://pulsar.princeton.edu/tempo) together with the JPL DE200 solar- 
system ephemeris ( |5tandish 1990|) to fit the measured pulse arrival times to the timing 
model with a least-squares technique. Results for the astrometric, spin, and dispersion 
parameters of PSR B1534+12 are presented in Table 2. Initial tests indicated the likely 
presence of systematic errors in the Arecibo data taken at 430 MHz (see below), so these 
data were used only for the calculation of the dispersion correction. The uncertainty quoted 
for DM is dominated by the unavoidable difficulty in aligning the frequency-dependent pulse 
shapes obtained in the different observing bands. We list for the first time a significant 
measurement of the rate of change of dispersion measure for this pulsar, determined from 
the dual-frequency data collected at Arecibo from 1990 to 1994. We note that there is no 
guarantee that this parameter will have remained constant since 1994. 

We fit the data to two models of the pulsar orbit. The theory-independent "DD" model 
( Pamour fc Deruelle 19861 ) treats all five PK parameters defined in §3.1 as free parameters in 



the fit. Alternatively, the "DDGR" model ( [laylor 19871 , [lay lor fc Weisberg 1989J) assumes 



general relativity to be correct and uses equations |] through [12] to link the PK parameters 
to M = mi +m 2 and m^-, consequently it requires only two post-Keplerian free parameters. 

Table |3| presents our adopted orbital parameters. Uncertainties given in the table are 
approximately twice the formal "1 a" errors from the fit; we believe them to be conservative 
estimates of the true 68%-confidence uncertainties, including both random and systematic 
effects. The Keplerian orbital parameters include the period Pb, projected semi-major axis 
x = aisini/c, eccentricity e, longitude of periastron u, and time of periastron To. These 
quantities are followed by the measured post-Keplerian parameters relevant to each of the 
two models. For the first time, the precision of this experiment requires the DDGR solution 
to include a parameter we call "excess P&," which accounts for an otherwise unmodeled 
acceleration resulting from galactic kinematics. 

The best estimates of the masses of the pulsar and its companion come from the DDGR 
solution. We find the masses to be equal: m x = m 2 = 1.339 ± 0.003 M Q . For the sake of 
comparison Table |3] lists, in italic numbers, computed PK parameter values corresponding 
to the measured masses in the DDGR fit. We call attention to the fact that the fitted and 
derived parameter values are in excellent accord, indicating good agreement of the theory- 
independent solution with general relativity. 

Figure 1 shows the post-fit residuals for the Arecibo 1400 MHz data, the Green Bank 
data, and the Jodrell Bank data, plotted as functions of date. Even within a single data set, 
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Fig. 1. — Post-fit residuals versus date for (a) Arecibo f 400 MHz, (b) Green Bank and (c) 
Jodrell Bank data. 

the TOA uncertainties can vary by factors of three or more, above and below the median 
values o"toa listed in Table 1, because of scintillation-induced intensity variations. We have 
not attempted to show these differences in data quality by means of error bars in the residual 
plots. Figure 2 illustrates the average post-fit residuals for the Arecibo measurements at 
1400 MHz and the Jodrell Bank data at 610 MHz, plotted as functions of orbital phase. 



3.3. Data Confidence Tests 



For some time we have suspected that solutions for the fitted parameters of PSR 
B1534+12, as determined from the Arecibo 430 MHz data, are biased by small systematic er- 
rors in the TOAs. These errors likely arise from imperfect post-detection dispersion removal 
that follows as a consequence of the relatively coarse frequency resolution of the filter-bank 
spectrometer. It is extremely difficult to achieve TOAs with systematic errors less than a 
few percent of the dominant instrumental smoothing effects in the data — in this instance, 
the dispersion time across a single channel. Variable spectral features caused by interstellar 
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Fig. 2. — Average post-fit residuals as a function of orbital phase for (a) the Arecibo 
1400 MHz data and (b) the Jodrell Bank 610 MHz data. 



scintillation, together with slightly irregular filter passbands and center frequencies, make 
the Arecibo 430 MHz measurements unreliable at the several microsecond level. 

Some further information on the data quality is presented in Figure 3. Here we follow 
Taylor and Weisberg (1989) and explore the statistical properties of the post-fit residuals to 
see if they "integrate down" as "nT 1 ! 2 when n consecutive values are averaged. As expected, 
this test confirms that the Arecibo 430 MHz data (and to a lesser extent the Green Bank 
data, which anyway receive low weight) deviate from the ideal rr x l 2 behavior. This indicates 
the likely presence of systematic errors in the TOAs, though we cannot completely rule out 
the possibility of a frequency- dependent deviation from the timing model (e.g., aberration- 
induced pulse profile variations). In any case, it is clear that for the highest accuracy within 
our preferred models we must omit the existing Arecibo 430 MHz measurements from our 
solution, as described above. 

The parameters of PSR B1534+12 are therefore best determined from the Arecibo 
1400 MHz data along with the Green Bank and Jodrell Bank data sets. The higher ob- 
serving frequency of 1400 MHz assures that for these measurements the total instrumental 
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Fig. 3. — Root-mean-square residual versus number of consecutive residuals averaged, for 
the four data sets: Green Bank, Jodrell Bank, and Arecibo at 430 and 1400 MHz. For the 
Arecibo 430 MHz data, filled circles represent the DD solution presented in Table 3, and 
open circles represent a DDGR fit to the 430 MHz data alone. Dashed lines indicate the 
expected slope of —1/2 for uncorrelated residuals. 



smoothing, even with 1.25 MHz channel bandwidths, is only 97/xs (see Table 1). Better still, 
the Princeton Mark IV system used at Jodrell Bank was explicitly designed to minimize or 
eliminate dispersion-related systematic errors in TOAs. We used these two data sets, both 
with and without the much lower-weighted Green Bank data, to measure and check the sys- 
tem parameters with highest accuracy. Arbitrary offsets were allowed between the different 
data sets to allow for frequency-dependent changes in pulse shape and slight differences in 
the standard profile alignments. For obvious reasons, we look forward to making improved 
measurements at 430 MHz with the Princeton Mark IV system, when the Arecibo telescope 
is available for use following its upgrade. 

Note that each of the data sets treated in Figure 3 may deviate at least slightly from 
the ideal n -1 / 2 behavior, suggesting small systematic errors of unknown origin. To partially 
compensate for this effect, uncertainties in TOAs were increased (in quadrature) by 2.9, 17.1, 
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and 20.5 /is, respectively, when calculating weights for the Arecibo, Jodrell Bank, and Green 
Bank data sets. 



4. Discussion 

4.1. Observed Change in Orbital Period: Distance to PSRB1534+12 

Before proceeding to tests of GR we must correct the observed P&, which is expressed in 
the reference frame of the solar system barycenter, to the center-of-mass frame of the binary 
pulsar. Damour and Taylor (1991) derived an expression for the most significant bias, which 
arises from galactic kinematic effects. It can be written as the sum of terms arising from 
acceleration toward the plane of the Galaxy, acceleration within the plane of the Galaxy, 
and an apparent acceleration due to the proper motion of the binary system: 




gal 



a 7 sin b 



cos I H — — cos A 



cRq cR\ 



+ A (13) 



Here a z is the vertical component of galactic acceleration, I and h the galactic coordinates of 
the pulsar, /i the proper motion, w the circular velocity at the Sun's galactocentric radius 
Ro, vi and R\ the corresponding values at the pulsar's location, A the angle between the 
Sun and the galactic center as seen from the pulsar, and d the distance from the pulsar to 
the Sun. The pulsar distance can be estimated from the dispersion measure, together with a 
smoothed-out model of the free electron distribution in the Galaxy. The Taylor and Cordes 
(1993) model yields d ~ 0.7kpc for PSR B1534+12, with an uncertainty of perhaps 0.2kpc. 
At this distance we estimate a z /c = (1.60 ± 0.13) x 10 -19 s _1 from the model of Kuijken and 
Gilmore (1989). Following Damour and Taylor (1991), we assume a flat galactic rotation 
curve and take vq = vi = 222±20kms _1 and Ro = 7.7±0.7kpc. Then, summing the terms 
in equation ( |i~3"|) and multiplying by P&, we find the total kinematic correction to be 

(iV) 8al = (0.038 ± 0.012) x 10~ 12 . (14) 

The uncertainty in this correction is dominated by the uncertainty in distance, which is only 
roughly estimated by the Taylor and Cordes model. 



Our measurement of the intrinsic rate of orbital period decay is therefore 

'p fe ) gal = (-0.167 ± 0.018) x 10~ 12 . (15) 



(A) 



Under general relativity, the orbital period decay due to gravitational radiation damping, 
(P,) GR , can be predicted from the masses mi and rri2 (eq. [10]), which in turn can be deduced 
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from the high precision measurements of u and 7. The expected value is 



(A) 



GR = -0.192 x KT 12 . (16) 



Although the measured value in the pulsar center-of-mass frame differs from this prediction 
by 1.4 standard deviations, it can be brought into good agreement by increasing the pulsar 
distance to slightly over lkpc. Stated another way, we can assume that GR is the correct 
theory of gravity, measure the "excess P&" for the system as described above and presented 
in Table 3, and then invert equation O to determine the pulsar distance (|Bell fc Bailes 



1996| ). By this method we obtain d = 1.1 ±0.2 kpc (68% confidence limit). The uncertainty 



is dominated by the measurement uncertainty of (P;,) obs , rather than uncertainties in the 
galactic rotation parameters or the acceleration a z . It follows that continued timing of this 
system should lead to a much more precise distance measurement, in due time. 

We note that the timing solution provides a second, independent constraint on the 
distance. The upper limit on parallax, it < 1.7mas (Table 2) constrains the distance to 
d > 0.6 kpc. While the parallax distance has less precision than the kinematic distance, it is 
reassuring that these measurements are in agreement. 

An accurate distance for PSR B 1534+ 12 is of considerable interest because of the im- 
portance of this system to estimates of the rate of coalescence of binary neutron-star pairs 
in a typical galaxy. Previous estimates of this rate have used much smaller distances for 
PSR B1534+12, including 0.4 kpc fllNarayan, Piran, fc Shemi 1991Q , 0.5 kpc QPhinney 19911 ), 



and 0.7 kpc flCurran fc Lorimer 1995| , [van den Heuvel fc Lorimer 1996| ), leading to a low 



estimate of its intrinsic luminosity. If our newly determined distance is correct, as we be- 
lieve it must be, the earlier studies have overestimated the number density of systems like 
PSR B1534+12 in the universe by factors of 2.5 to 20. We find that the most recent obser- 
vational estimate of the double neutron star inspiral rate, 2.7 x 10~ 7 yr _1 ([van den Heuvel 15 



Lorimer 19961 ), must be reduced to between 1.1 x 10 7 y r X an d 1.4 x 10 7 yr , depending on 



the unknown scale height of such systems in the Galaxy. This is, of course, a lower bound 
on the number of coalescence events; the actual rate will be higher by unknown factors that 
account for the fraction of pulsars with beams that do not cross the Earth and the number 
of double neutron star binaries that do not contain an active pulsar. 



4.2. Test of Relativity 

Our analysis of this experiment represents the second test of general relativity based on 
the ui, 7, and Pb parameters of a binary pulsar system, and the first one to add significant 
measurements of the Shapiro-delay parameters r and s. The left-hand sides of equations 
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Fig. 4. — Mass-mass diagram for the PSR B1534+12 system. Labeled curves illustrate 
68% confidence ranges of the DD parameters listed in Table 3. The filled circle denotes the 
component masses according to the DDGR solution. A kinematic correction for assumed 
distance d = 0.7 ± 0.2 kpc has been subtracted from the observed value of P^. A slightly 
larger distance removes the small apparent discrepancy. 



12]) represent measured quantities, as specified for this experiment in the "DD" column of 
Table 3. If GR is consistent with the measurements and there are no significant unmodeled 
effects, we should expect the five curves corresponding to equations (PHi"2"D to intersect at a 
single point in the m 1 -r7i 2 plane. A graphical summary of the situation for PSR B1534+12 
is presented in Figure [|, in which a pair of lines delimit the 68% confidence limit for each 
PK parameter (a single line for u, whose uncertainty is too small to show). A filled circle 
at rrii — m 2 — 1-339 M Q marks the DDGR solution of Table 3, and its location on the u 
line agrees well (better than 1%) with the measured DD values of 7 and s. We have already 
noted that the DD value of Pb is slightly too small when corrected to the dispersion-estimated 
distance. However, as discussed above, this discrepancy can be removed by invoking a larger 
distance to the pulsar. Finally, the value of r, although presently little better than a 20% 
measurement, is also well centered on its expected value. Altogether we find the measured 
parameters of the PSR B1534+12 system to be in excellent accord with general relativity, 
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and consequently this theory has passed another very significant astrophysical test. 
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Table 1: Parameters of the four observing systems. 



Arecibo Arecibo Green Bank Jodrell Bank 

Mark III Mark III Spec. Proc. Mark IV 



Frequency (MHz) 430 1400 

Bandwidth (MHz) 8 40 

Spectral Channels 32 32 

Dedispersing system . . . incoherent incoherent 

Time resolution (/xs) . . . 329 97 

Integration time (s) ... 180 300 

Dates 1990.7-94.2 1990.8-94.1 

Number of TOAs 2311 1170 

Median ct TO a (^ s ) 4.9 5.4 



575, 800 


610 


40 


5 


512 


1 


incoherent 


coherent 


296 


19 


180 


190 


1994.2-97.4 
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685 


780 


40 
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Table 2: Astrometric, spin, and dispersion parameters for PSR B1534+12? 

Right ascension, a (J2000) 15 h 37 m 09!95994(2) 

Declination, 5 ( J2000) 1 1° 55' 55'/6561 (3) 

Proper motion in R.A., \x a (masyr -1 ) .... 1.3(3) 

Proper motion in Dec, [x& (masyr -1 ) —25.5(4) 

Parallax, n (mas) < 1.7 

Pulse period, P (ms) 37.9044404878552(5) 

Period derivative, P (10~ 18 ) 2.42253(3) 

Epoch (MJD) 48778.0 

Dispersion measure, DM (cm~ 3 pc) 11.619(12) 

DM derivative (cm-^cyr" 1 ) -0.00036(2) 

Galactic longitude I (deg) 20.0 

Galactic latitude b (deg) \1.8 

Composite proper motion, \x (masyr -1 ) . . . 25.5(4) 

Galactic position angle of \x (deg) 238. 7(2) 



a Figures in parentheses are uncertainties in the last digits quoted, and italic numbers represent derived 
quantities. 
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Table 3: Orbital parameters of PSR B1534+12 in the DD and DDGR models a 



DD model 



DDGR model 



Orbital period, P b (d) 0.42073729929(4) 0.42073729930(4) 

Projected semi-major axis, x (s) 3.729463(3) 3.7294628(7) 

Eccentricity, e 0.2736777(5) 0.2736776(2) 

Longitude of periastron, u (deg) 267.44746(16) 267.44760(12) 

Epoch of periastron, T (MJD) 48777.82595097(6) 48777.82595096(6) 

Advance of periastron, Cj (degyr -1 ) 1.75576(4) 1.75577 

Gravitational redshift, 7 (ms) 2.066(10) 2. 057 

Orbital period derivative, (P 6 ) obs (10~ 12 ) . . . -0.129(14) -0.1924 

Shape of Shapiro delay, s 0.982(7) 0.9797 

Range of Shapiro delay, r (yus) 6.7(1.3) 6.60 

Derivative of x, \x\ (10~ 12 ) < 2.5 < 0.022 

Derivative of e, |e| (10 -15 s _1 ) <6 <6 

Total mass, M = m x + m 2 (M ) 2.67838(8) 

Companion mass, m 2 (M ) 1.339(3) 

Excess P b (10~ 12 ) 0.062(14) 



a Figures in parentheses are uncertainties in the last digits quoted. Italic numbers represent derived parameters, 
assuming general relativity. 



